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Abstract 

The dynamical equations of clarinet-like systems are known to be re- 
ducible to a non-linear iterated map within reasonable approximations. 
This leads to time oscillations that are represented by square signals, 
analogous to the Raman regime for string instruments. In this article, we 
study in more detail the properties of the corresponding non-linear itera- 
tions, with emphasis on the geometrical constructions that can be used to 
classify the various solutions (for instance with or without reed beating) 
as well as on the periodicity windows that occur within the chaotic region. 
In particular, we find a regime where period tripling occurs and examine 
the conditions for intermittency. We also show that, while the direct ob- 
servation of the iteration function does not reveal much on the oscillation 
regime of the instrument, the graph of the high order iterates directly 
gives visible information on the oscillation regime (characterization of the 
number of period doubligs, chaotic behaviour, etc.). 

Keywords ; Bifurcations, Iterated maps. Reed musical instruments, Clarinet, 
Acoustics. 



1 Introduction 

Non-linear iterated maps are now known as an universal tool in numerous sci- 
entific domains, including for instance mechanics, hydrodynamics and econ- 
omy [I] [5] [3]. They often appear because the differential equations describing 
the dynamics of a system can be reduced to non-linear iterations, with the 
help of Poincare recurrence maps for instance. The resulting iterations combine 
a great mathematical simplicity, which makes them convenient for numerical 
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simulations, with a large variety of interesting behaviors, providing generic in- 
formation on the properties of the system. In particular, they are essential to 
characterize one of the routes to chaos, the cascade of period doublings [1]. 

In musical acoustics, Mc Intyre et al. have given, in a celebrated article [S], a 
general frame for calculating the oscillations of musical instruments, based upon 
the coupling of a linear resonator and a non-linear excitator (for reed instru- 
ments, the flow generated by a supply pressure in the mouth and modulated by 
a reed). In an appendix of their article they show that, within simplified models 
of self-sustained instruments, the equations of evolution can also be reduced to 
an iterated map with appropriate non-linear functions. For resonators with a 
simple shape such as a uniform string or a cylindrical tube, the basic idea is to 
choose variables that are amplitudes of the incoming and outgoing waves (trav- 
elling waves) , instead of usual acoustic pressure and volume velocity in the case 
of reed instruments. If the inertia of the reed is ignored (a good approximation 
in many cases), and if the losses in the resonator are independent of frequency, 
the model leads to simple iterations; the normal oscillations correspond to the 
so called "Helmholtz motion" , a regime in which the various physical quantities 
vary in time by steps, as in square signals. Square signals obviously are a poor 
approximation of actual musical signals, but this approach is sufficient when the 
main purpose is to study regimes of oscillation, not tone-color. 

In the case of clarinet- like systems, the idea was then expanded 0, giving rise 
to experimental observations of period doubling scenarios and to considerations 
on the relations between stability of the regimes and the properties of the second 
iterate of the non-linear function; see also [7] and especially [5] for a review 
of the properties of iterations in clarinet-like systems and a discussion of the 
various regimes (see also [9]). More recent work includes the study of oscillation 
regimes obtained in experiments |10llll| . computer simulation jl2| as well as 
theory [T1[T1]. 

The general form of the iteration function that is relevant for reed musical 
instruments is presented in section |31 It it is significantly different from the 
usual iteration parabola (i.e. the so-called logistic map). Moreover, it will be 
discussed in more detail that the control parameters act in a rather specific way, 
translating the curve along an axis at 45° rather than acting as an adjustable 
gain. 

The purpose of the present article is to study the iterative properties of 
functions having this type of behavior, and their effect on the oscillation regimes 
of reed musical instruments. We will study the specificities and the role of the 
higher order iterates of this class of functions, in particular in the regions of 
the so called "periodicity windows" , which take place beyond the threshold of 
chaos. These windows are known to contain interesting phenomena f^|15516j. for 
instance period tripling or a route to intermittence, which to our knowledge have 
not yet been studied in the context of reed musical instruments. Moreover, the 
iterates give a direct representation of the zones of stability of the different 
regimes (period doublings for instance), directly visible on the slope of the 
corresponding iterate. 

For numerical calculations, it is necessary to select a particular representa- 
tion of the non-linear function, which in turn requires to choose a mathematical 
expression of the function giving the volume flow rate as a function of the pres- 
sure difference across the reed. A simple and realistic model of the quasi-static 
flow rate entering a clarinet mouthpiece was proposed in 1974 by Wilson and 
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Beavers [T7], and discussed in more detail in 1990 by Hirschberg et al. [H]. This 
model provides a good agreement with experiments |19j and leads to realistic 
predictions concerning the oscillations of a clarinet |20] . Using this mathemati- 
cal representation of the flow rate, we will see that iterations lead to a variety of 
interesting phenomena. Our purpose here is not to propose the most elaborate 
possible model of the clarinet, including all physical effects that may occur in 
real instruments. It is rather to present general ideas and mathematical solutions 
as illustration of the various class of phenomena that can take place, within the 
simplest possible formalism; in a second step, one can always take this simple 
model as a starting point, to which perturbative corrections are subsequently 
added in order to include more specific details. 

We first introduce the model in § [21 and then discuss the properties of the 
iteration function in § [31 The bifurcations curves are obtained in § [H and, in § 
[31 we discuss the iterated functions and their applications in terms of period 
tripling and intermittence. In particular we see how the graph of high order 
iterates give visible information on the regime of oscillation (number of period 
doublings for instance) or the appearance of a chaotic regime, while nothing 
special appears directly in the graph of the first iterate. Two appendices are 
added at the end. 



We briefly recall the basic elements of the model, the non-linear characteristics 
of the excitator, and the origin of the iterations within a simplified treatment 
of the resonator. 

2.1 Nonlinear characteristics of the entering flow 

In a quasi static regime, the flow U entering the resonant cavity is modelled 
with the help of an approximation of the Bernoulli equation, as discussed e.g. 
in [15]. We note Pint the acoustic pressure inside the mouthpiece, assumed to 
be equal to the one at the output of the reed channel, P,„ the pressure inside 
the mouth of the player; for small values of the difference: 



the reed remains close to its equilibrium position, and the conservation of en- 
ergy implies that U is proportional to 'r]p^y\AP\, where 77^ = ±1 is the sign 
of AP (we ignore dissipative effects at the scale of the flow across the reed 
channel) ; for larger values of this difference, the reed moves and, when the dif- 
ference reaches the closure pressure Pc, it completely blocks the flow. These two 
effects are included by assuming that if AP < Pc the flow U is proportional 
to rjp^yiAPl [Pc — AP], and if AP > Pc, the flow vanishes. Introducing the 
dimensionless quantities: 



2 The model 




(1) 



P = Pint/Pc 
U = UZoo/Pc 
7 - P^/Pc 



(2) 



Ap = AP/Pc = 7 - P • 
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where Z^q = pcjS is the acoustic impedance of an infinitely long cylindrical 
resonator having the same cross section S than the clarinet bore (p is the density 
of air, c the velocity of sound), we obtain: 

w = (3) 
if Ap > 1 i.e. p < 7 — 1 ; 

u = C{P + ^ - - P (4) 

if < Ap < 1 i.e. 7 - 1 < p < 7 ; 

u = -({p+l - 7)\/p- 7 (5) 
if Ap < i.e. p > "f. 

The parameter ^ characterizes the intensity of the flow and is defined as: 

where Sop is the opening cross section of the reed channel at rest. One can show 
that C is inversely proportional to square root of the reed stiffnesf0, contained 
in Pc- In real instruments, typical values of the parameters are 7 € [0,1.5]; 
C = [0.1, 0.5] ; values C > 1 will not be considered here, since they correspond to 
multi- valued functions u{p), a case that does not seem very realistic in practice. 
Fig IT] shows an example of function defined in Eqs.(l3] to [5]). It is obviously 
non-analytic; it is made of three separate analytic pieces, with a singular point 
at p = 7. The derivative of the function u{p) is discontinuous at p = 7 — 1 
(point Mb in Fig. [U the index b being used for the limit of possible beating); 
a smoothing of the resulting angle of the function could easily be introduced at 
the price of a moderate mathematical complication, but this is not necessary 
for the present discussion. 



2.2 Iteration 

Waves are assumed to be planar in the quasi one dimensional cylindrical res- 
onator. Any wave can be expanded into an outgoing wave p'^(t — z/c) and an 
incoming wave p^ (t + z/c), where t is the time and z the abscissa coordinate 
along the axis of the resonator; at point z = (at the tip of the reed), the 
acoustic pressure and flo'nQ are given by: 

p{t)^p+(t)+p-it) ; u{t)^p+{t)-p-{t) (7) 



P+{t)^-[p{t)+u{t)] ; p-{t)^-[p{t)~u{t)]. (8) 

We will use variables p^(i) instead of p{t) and u{t). If we assume that the 
impedance at the output of the resonator is zero (no external radiation, the 

^the reed remains close to its equilibrium position; the acoustic flow is then independent 
of the stiffness of the reed. Equation ||4} then provides UZoo/Pc — C \/ {Pm — Pint)/ Pc, or 
U ~ {CVPc/Zcxd) \J {Pm — Pint)', but Pc is roughly proportional to the reed stiffness, so 
that the independence of the flow with respect to the stiffness requires that is inversely 
proportional to the square root of this stiffness. 

^The flow is related to the pressure via the Euler equation: dP/dz = —pS^^dU/dt. 
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Figure 1: Graph showing the air flow rate u entering the resonator of a clarinet as 
a function of the internal pressure p (pressure in the mouthpiece) . All physical 
quantities are expressed in dimensionless units, as explained in the text. Mf, 
corresponds to the contact point where the internal pressure in the mouthpiece 
bends the reed sufficiently to close the channel, so that the flow vanishes; it 
remains zero in all region p < 7 — 1. is the inversion point where p = 7 and 
where the acoustic flow changes sign. The full line corresponds to 7 0.4, the 
broken line to 7 = 0.6. Here C, = 0.8. 



output pressure remains the atmospheric pressure), we obtain the reflection 
condition: 

p-(t) = -p+(t-2£/c) , (9) 

where i is the resonator length and c the sound velocity. This equation expresses 
that the reflected wave has the same amplitude than the incoming wave. Losses 
are not included in this relation, but one can also introduce them very easily by 
replacing ([9]) by: 

p-(t) = -Ap+(t-2£/c) , (10) 

which amounts to introducing frequency independent losses; a typical value is 
A — 0.9. For a cylindrical, open, tube with no radiation at the open end so that 
losses only occur inside the tube, A — exp(— 2a^), where a is the absorption 
coefficient. Of course this is an approximation: real losses are frequency depen- 
denlH and radiation occurs but, since losses remain a relatively small correction 
in musical instruments, using Eq. (jlO|) is sufficient for our purposes. 

We now assume that all acoustical variables vanish until time t = 0, and 
then that the excitation pressure in the mouth suddenly takes a new constant 
value 7; this corresponds to a Heaviside step function for the control parameter. 
Between time and time 2Z/c, according to (llOp . the incoming amplitude p^{t) 
remains zero, but the outgoing amplitude p'^{t) has to jump to value p'^ in order 

^The value of a depends on both frequency / and radius R. For normal ambient conditions 
(20°C), a = 2.96 IQ-'^^/R (see e.g. EI]). 
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to fulfil Eqs.(I3]to[S]) . At time t — 21/ c, the variable p~ [t) jumps to value — ApJ, 
which immediately makes p^{t) jump to a new value pj, in order to still fulfil 
Eqs.([3]to[5]). This remains true until time t = ^l/c, when p~{t) jumps to value 
—\P2 and to a value , etc. By recurrence, one obtains a regime where 

all physical quantities remain constant in time intervals Inljc < t < 2(n + l)l/c, 
in particular p„ for the pressure and w„ for the flow, with the recurrence relation: 

Pn - (11) 

In what follows, it will be convenient to use 21/c as a natural time unit. We will 
then simply call "time n" the time interval {n—l)2l/c<t< n2l/c. Notice that 
in order to get higher regimes (with e.g. triple frequency), the previous choice 
of transient for 7 needs to be modified (see e.g. [5]). 

Now, by combining Eqs.(|3]to[5l) and [7]), one can obtain a non- linear relation 
g between p^ and p^ : 

vl^giPn), (12) 

which, combined with (jlip . provides the relation: 

P;l-.9(-Ap+_i) = /(p+_i), (13) 

with, by definition: f{x) = g{—\x). The equation of evolution of the system are 
then equivalent to a simple mapping problem with an iteration function f{x). 
The graph of this function is obtained by rotating the non-linear characteristics 
of Fig. [Uby 45° (in order to obtain g), then applying a symmetry (to include the 
change of sign of the variable) and finally a horizontal rescaling by a factor 1/A; 
the result is shown in Fig. [2j This provides a direct and convenient graphical 
construction of the evolution of the system [B]; Fig. [3]shows how a characteristic 
point 1 is transformed into its next iterate 2, etc... by the usual construction, 
at the intersection of a straight line with the iteration curve, i.e. by transferring 
the value of f{x) to the x axis and reading the value of the function at this 
abscissa in order to obtain / [f{x)]. 

In what follows, we consider 7 as the main control parameter of the iteration; 
it corresponds to a change of pressure in the mouth of the instrumentalist. A 
second control parameter is which the player can also change in real time 
by controlling the lip pressure on the reed. For a given note of the instrument, 
parameter A remains fixed, but of course depends on which lateral holes of the 
clarinet are closed, in other words on the pitch of the note. 

The oscillations where the functions remain constant and jump to a different 
value at regular interval of times are reminiscent of the Raman regime for the 
oscillation of bowed strings [32]. Mc Intyre et al. have indeed noticed that, if 
one replaces the non-linear function by that corresponding to a bowed string, 
one obtains the Raman oscillation regime of a string bowed at its center [5]. 

3 Properties of the iteration 

The analytical expression of the iteration function is given in Appendix [X] 
Figure [5] shows the function for given values of the parameters 7 and C, and 
three different values of the loss parameter A. 

In the literature, the most commonly studied functions have the following 
properties (see e.g. Collet and Eckman |3I or Berge et al [2]): 
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Figure 2: Iteration function / for C — 0.8 and 7 = 0.43 Solid line A = 1 (no loss 
in the resonator); mixed line A = 0.8; dotted line A = 0.4. The circles on the 
right indicate the contact point Mf,, those on the left the flow inversion point 



• They are defined on a finite interval and map this interval into itself; 

• They are continuous; 

• They have a unique maximum; 

• their Schwarzian derivative is negative. 

A function verifying these properties will be called a "standard" function; 
the function f{x) of interest in our case does not fulfil all these requirements. 

Domain of iteration Usually, the iteration function defines an applica- 
tion of the interval < x < 1 over itself. Here, /(x) is defined on an infinite 
interval [—00, +00] even if, obviously, very large values of the variables are not 
physically plausible. Nevertheless, analyzing the different cases corresponding 
to Eqs.(l3] to [5]), one can show that the function f{x) has a maximum /max 
obtained for: 




(14) 



with value: 



max 



7: + Aq with 




) 



18 



1 



(15) 



where x is defined by: 




(16) 



It can be shown that this maximum is unique for large value of C (C > 1 /•s/S); 
for smaller values, a second maximum exists at a very large negative values of 



i.e. for very large negative flow, but we will see below that such values of the 
flow cannot be obtained after a few iterations. Therefore we focus our attention 
only on the maximum /max, which varies slowly as a function of C because 
increases monotonically from for C = to a small value (5/54, for C = 1). 

The geometrical construction of Fig. [3] shows that, after a single iteration, 
the characteristic point M necessarily falls at an abscissa x < /max- Let us call 
/min = fifmax) the Ordinate of the point on the iteration function with abscissa 
/max- The two vertical lines x = /min and x = /max, together with the two 
horizontal lines y = /min and y = /max, define a square in the x, y plane, from 
which an iteration cannot escape as soon as the iteration point has fallen inside 
itB Conversely, since every characteristic point has at least two antecedents, the 
iteration can bring a point that was outside the square to inside. In other words, 
the square determines a part of the curve which is invariant by action of the 
function. For usual initial conditions, such as Pq = 0, the starting point already 
lies within the square, so that all points of the iteration keep this property. We 
have checked that, even if one starts with very large and unphysical pressure 
differences (positive or negative), the iterations rapidly converge to the inside 
the square. In what follows, we call it the "iteration square" . 

The net result is that, if we do not consider transients, we can consider that 
the function defines an application of the interval [/min, /max] over itself. We 
are then very close to the usual mapping situation, except that here the interval 
depends on the control parameters (since the value of /max depends on 7 and 
C), but with a relatively slow variation. 



Singularities An interesting feature of the iteration function is the discon- 
tinuity of its first derivative occurring at the beating limit point Mb at x — Xh, 
given by: 

= ; f(xb) = . (17) 

When the reed closes the channel (p+ — , p'^ + p^ < 7 — 1), x > Xb, 
f{x) ~ —Xx, the iteration function is linear. 

Another singularity, i.e. a discontinuity of the second derivative, is obtained 
at the crossover between positive and negative flow, the inversion point Mi 
where sign of the flow changes. Its abscissa x = Xi is given by: 

^. = -^ ; fi^d^l- (18) 

For < 7 < 1, is negative and Xb positive: therefore the initial point of 
the iteration (x — 0) lies in the interval [xi,Xf,], with neither contact with the 
mouthpiece nor negative flow, as one could expect physically. 



Schwarzian derivative The Schwarzian derivative 3 of f{x) is equal to: 

2 

(19) 



f_ 

f 



*We assume that f(fmin) > fmin, which means that the iteration curve crosses the left 
side of the square, as is the case in Fig[3l 
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Figure 3: Graphical illustration of the iteration, where an initial point 1 is 
iterated into point 2, 3, etc. (similarly for the point 1', 2', 3',.--)- Since the non- 
linear iteration function has a maximum fmax, after a few steps the iteration 
remains inside an "iteration square" shown in broken lines. This square has 
its upper side tangent to the maximum of the function, at point M^axi which 
after one iteration becomes point Mmin defining the lowest side of the square 
(ordinate fmin)- Depending on the parameters, the iteration square contains or 
does not contain the contact point Mf, and the flow inversion point Mi. Here 
7 = 0.44, C = 0.8, A = 0.95. 



where /', /" and /"' indicate the first, second and third derivatives of f{x), 
respectively. If a; > x;,, it is zero; ii Xi < x < Xb, using the change of variables 
given in Appendix 1X1 Sf can be shown to be equal to: 

o \2 

Sf - Y'\Y'-2)^ [Y"'Y\Y' - 2) - iY"\Y' - 1)] , (20) 

where Y is a function oi X - see Eqs. (|27l) to ([^5]) . Therefore its sign does 
not depend on the loss parameter A. After some calculations, the Schwarzian 
derivative is found to be negative for all x € [xi^ Xb] when C < 1/ y/E. Otherwise, 
it is negative up to a certain value, then positive \xg to x ~ xi,- The calculation 
of S f for the case x < Xi shows that it is positive, except for a small interval. 
The iteration function therefore differs from a standard function because of the 
sign of the Schwartzian derivative; this is related to the nature of the bifurcation 
at the threshold of oscillation [53], which can be either direct or inverse. 

Beating and negative flov^r limits In Fig. |3] we see that, depending 
whether the contact point Mb and flow inversion point Mi of the iteration curve 
fall inside or outside the iteration square, a beating behavior of the reed and a 
sign inversion of the air flow are possible or not. 

Point Mb falls inside the iteration square if its abscissa Xb given by ([T7)l is 
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smaller than fmax, which leads to: 

7 > 7. ^ ^TX^- (21) 
The limiting value is less than unity (it tends to 1/2 when tends to and 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 4: In the plane of the control parameters 7 and ^, the line representing 
7h gives a border between the upper region, where reed beating may occur, and 
the lower region where it cannot - see Eq. (I2ip . As a point of comparison, the 
line labelled 7f,2s correspond to the limit obtained in [14] for the particular case 
of a 2-state regime, and given by Eg. d^^ . The figure also shows the line and 
7i2 associated with the possibility of negative flow (the flrst one turns out to 
be very close to that associated with reed beating). Small values for the losses 
have been assumed (A — 0.95). 

A tends to unity). This necessary condition is completely independent of the 
nature of the limit cycle, and less stringent than the limit 7^25 obtained in [14) . 
for a 2-state cycle : 

762s = 1/2 [l + /3/3i+/3?(l-/3/3i)], (22) 

where /3 = C(l — A)/(l + A) and /3i = Fig. [5]gives a comparison between 

the two limits. Similarly, a necessary condition for possible inversions of the 
sign of the air flow is that point Mi falls inside the iteration square of Fig |31 
in other words that Xi is larger than fmin- We show in Appendix |B] that this 
happens if: 

7ii < 7 < 1^2■ (23) 

The expression of the two limits 7^1 and 7^2 are given in the Appendix and 
can be seen on Fig. [H They are solutions of Xi — fmim and exist only if the 
following condition holds: 
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Therefore, for a given C, negative flow is possible only above a certain value 
of A; this value is 27/32 = 0.84 for C = li tends to unity when ( tends to 
zero. Using a more realistic shape for the function f{x) with a rounding of the 
kink at Xb (no discontinuity of the derivative) should lead to a shorter range of 
negative flow, making the phenomenon even less likely, as illustrated in Fig. |4l 
Of course, the two above conditions (PT|) and (1^ are necessary, but not 
sufficient; they do not ensure that either beating or flow inversion will indeed 
take place, since this will be true only if the corresponding regions of the non- 
linear curve are reached during the iterations. Generally speaking, this will have 
more chance to occur in chaotic regimes, where many points are explored in the 
iterations, than in periodic regimes. Since no observation of negative flow has 
been reported in the literature, it is not clear whether this actually happens in 
real instruments. 

In conclusion of this section, the iteration function is similar to those usually 
considered in the context of iterated maps, without really belonging to the 
category of "standard" functions. The major difference is actually the effect of 
the control parameters on the function, since usually the control parameters 
acts as a gain, expanding the vertical axis of the graph; here the parameter 7 
(pressure in the mouth of the instrumentalist) translates the iteration function 
along an axis at 45° of the coordinate axis, while the other control parameter C 
(the pressure of the lip on the reed) expands the function along the perpendicular 
axis. It is therefore not surprising that we should find a parameter dependence of 
the dynamical behaviors that is significantly different from the standard results. 

4 Bifurcation curves 

Figure [5] shows an example of bifurcation curves, for A = 0.95 and ( = 0.8, 
and illustrates the relative complexity of the possible regimes. The upper curve 
corresponds to the outgoing amplitude (or x), the middle curve to acoustic 
pressure p, and the lower curve to the acoustic flow u. The three curves show 
the last 20 values obtained after computing 400 iterations for each value of 
the mouth pressure 7. By calculating 2000 iterations for a given value of the 
parameter 7, we have checked that the limit cycle is then reached. Obviously 
this method leads to stable regimes only. 

When the control parameter 7 increases, the beginning of these curves follows 
a classical scenario of successive period doublings, leading eventually to chaos; as 
expected, high values of the parameters C and A favour the existence of chaotic 
regimes, as well as beating reed or negative flow. When 7 continues to increase, 
another phenomenon takes place: chaos disappears and is replaced by a reverse 
scenario containing a series of frequency (instead of period) doublings. We call 
this phenomenon a "backwards cascade" (in order to distinguish it from the 
usual "inverse cascade", which takes place within periodicity windows inside 
chaos [2]); this backwards cascade is a consequence of the speciflcities of the 
effect of the control parameter on the iteration function in our model, and of 
the particular shape of the iteration function (for instance a straight line beyond 
the beating limit point). As a matter of fact, different kinds of cascades have 
been studied in the literature (see e.g. [53] and [IS], in particular Fig. 5). 
In Fig. [5] the variations of 7 correspond to a "crescendo": for a given value 
of 7, the initial value for the iteration, p^ , is chosen to be equal to the last 
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Figure 5: Bifurcation curve for A = 0.95 and C. — 0.8. For increasing values of 
the abscissa 7 (bfowing pressure), corresponding to a crescendo, the curve show 
the values after 400 iterations of the outcoming wave p"*" (top) the pressure p in 
the mouthpiece (middle), and the volume flow u (bottom). Above 7 ~ 0.45, the 
flow can be negative and the reed can beat. The top figure also shows /max(7) 
and /mill (7) (mixed lines) associated to the "iteration square". 

value p^QQ obtained with the previous value of 7. But we have also studied the 
"decrescendo" regime and observed that, in the chaotic regimes, the plotted 
points differ from the crescendo points; on the other hand, they remain the 
same in the periodic regimes, indicating a direct character of the bifurcations (no 
hysteresis). We have found an exception to this rule: between values 7 = 0.5 and 
7 = 0.53, 2-state and 4-state regimes coexist, indicating an inverse bifurcation. 
Another inverse bifurcation, between a 2-state regime and a static regime, occurs 
beyond the limit of the figure, the two regimes coexisting between 7=1 and 
7 = 6.3544; this is not shown here (the shape of the curve can be found in 
Ref. [2], see upper Fig. 4). 

The two limits of the function f{x), /max and /min, are also plotted in the 
upper figure ( p'^) showing that, as expected, the corresponding values remain 
inside the iteration square (§|31)- In the figure at bottom, the results for the flow 
u exhibit lower limits for negative flow and for beating, which are very close 
to the theoretical limits, respectively jn = 0.4454 and = 0.4503, and are 
located within the chaotic regime. Negative flow disappears at the bifurcation 
between the 4-state and the 2-state regime, 7 — 0.5262, a much lower value 
than the higher limit for negative flow 7^2 = 1.189. 

Table 1 shows the critical values of 7 corresponding to changes of regime. Up 
to the first chaotic regime (7 — 0.4409), the behavior follows the usual period 
doubling cascade scenario. Between 7 — 0.4467 and 7 = 0.4479 a "periodicity 
windows" [2] is obtained, with 6-state, then 1 2-state and 24-state regimes (but 
no 3-state regime). Above the value 7 = 0.4409 for which chaos starts, an 
"inverse cascade" type scenario is observed, then intermittences occur, chaos 
again, and finally the "backwards cascade" to the static regime. We did not try 
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From 7 ~ 


Regime 


Comments 


From 7 — 


Regime 


Comments 





1-state 




0.4540 


24-state 




0.3545 


2-state 




0.4542 


12-state 




0.4272 


4-state 




0.4544 


6-state 




0.4384 


8-state 




0.466216 


I 




0.4403 


16-state 




0.4664 


chaos 




0.4408 


32-state 




0.46945 


60-state 




0.4409 


chaos 


J^l = 0.4454 


0.4695 


20-state 




0.4467 


6-state 


PW 


0.4696 


chaos 




0.4474 


12-state 


PW 


0.46985 


4-state 




0.4476 


24-state 


PW 


0.5000 


2-state (D) 




0.4479 


chaos 


7b = 0.4503 


0.53 


2-state (C) 




0.4538 


36-state 




1. 


1-state (D) 


1^2 = 1.189 


0.4539 


chaos 




6.3544 


1-state (C) 





Table 1: Values of the parameter 7 at the lower limit of the different regimes, 
corresponding to Fig. [5] I=intermittencies; C= crescendo; D= decrescendo; 
PW= periodicity windows. 



to obtain the same accuracy for the values of all different thresholds, because the 
ranges for 7 have very different widths; for some values of 7, it has been necessary 
to make up to 2000 iterations, and sometimes it is not obvious to distinguish 
between a chaotic regime, a long transient, or an intermittency regime. 

5 Iterated functions 

We now discuss how the iterated functions can be used to study the different 
regimes and their stability. We write f'^'^\x) the second iterate of /, and more 
generally /("^(x) its iterate of order n; the derivative of / with respect to x is 
fix). Around the fixed point x* of the first iterate f{x), a Taylor expansion 
gives: 

fix) = fix*) + ix- x*)f'ix*) + .. = x* + ix- x*)f'ix*) + .. , 
which provides the well known stability condition for a fixed point x* of fix): 

\f'ix*)\ < 1. (25) 
Since the derivative of the iterate of order n is given by: 

one can show by recurrence that, when a; = a;*, it is equal to the n-th power of 
the derivative of /, so that: 

f^-\x)^x* + ix-x*)[f'ix*)r + .. (26) 

If the fixed point is stable (resp. unstable) with respect to fix), it is also 
stable (resp. unstable) with respect to any iterate. If a; is a vector, instead of 
a scalar, this linearized approach leads to the Floquet matrix, and fix) should 
be replaced by the eigenvalues of the matrix. 
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5.1 Stability of the period doubling regimes 




Figure 6: Iteration functions for A ~ 0.95, 7 = 0.31 and C ~ 0.8. The 1st 
iterate /(x) is shown with a mixed hue, the 2nd iterate /'•^''(x) with a (blue) 
soHd Hue, and (x) with a (red), thin mixed Hue. The dotted hnes are the 
first diagonal and the straight line perpendicular at the fixed point M* of /(x), 
solution of /(x) = X. The tangent lines to iterate 1 at the point M* is shown 
with a solid line. 

Examples of iterated functions of order 1, 2, and 4 are shown in figures [H] 
and [71 with the same values of C and A as in figure [SJ in the former, the blowing 
pressure 7 is 0.31, in the latter, 7 is 0.42. The first iterate has a unique fixed 
point, M* = (a;*, x*), located by definition on the first diagonal. The fixed point 
is stable if the absolute value of the derivative at M* is smaller than unity, in 
other words if the tangent line lies between the first diagonal (with slope +1) 
and its perpendicular (with slope —1). When 7 — 0.31, we see in Fig.[S]that the 
fixed point M* is stable, so that no oscillation takes place. When 7 increases, 
M* becomes instable and, at the same time, gives rise to three fixed points of 
j(2)^ For 7 — 0.42, Fig. [7] shows that the tangent is outside the angle between 
the diagonal and its perpendicular, so that the fixed point is now unstable; on 

the other hand, the second iterate /f^-* now has two more fixed points M^^'' 
(2) 

and M2 with slopes less than 1 (in absolute value) : we therefore have a stable 
2-state regime. 

The same scenario then repeats itself when 7 continues to increase: at some 

(2) (2) 

value, points AI^ and become instable in turn (the corresponding slope 

exceeds 1 in absolute value), and both points M^"^^ and M2^^ divide themselves 
into three fixed points of /*^^^; the two extreme new points have small slopes 
for this iterate, which leads to a 4-state stable regime. By the same process of 
successive division of fixed points of higher and higher iterates, one obtains an 
infinite number of period doublings, until eventually chaos is reached. This is 
the classical Feigenbaum route to chaos. 
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Figure 7: Iterates for A — 0.95, 7 — 0.42 and C = 0.8, with the same plots that 
in figureini The tangent hnes at the new fixed points of /^^^(x), and M'^p , 
are also shown. 

Some general remarks are useful to understand the shape of the iterates in 
the figures: 

• If the value of /(cc) for the abscissa x verifies j{x) = f{x*), i.e. if the 
point M{x, f{x)) is on a horizontal line y = x* , all iterates go through 
the same point; 

• The extrema of f^^\x) verify either f'{x) — (i.e. x — x^a^) or /'■^•'(a;) = 
/max , because df^'^\x)/dx — f [f{x)] f'{x); therefore the extrema of 
f'^'^\x) are at either the same abscissa or the same ordinate as those of 

/(^); 

• More generally, for n > 1, if /(""^^(x) = Xmax, then f^'^\x) = /,„ax, and 
it is at a maximum (its first derivative vanishes and the second one is 
negative), and if /(""^^(x) = /max, then f^'^^{x) = /min, and it is at a 
minimum (its first derivative vanishes and the second one is positive); 

• The kink of the first iterate (beating limit point) is also visible on the 
iterates; 

• A well known property of the Schwarzian derivative is as follows : If the 
Schwarzian derivative of f{x) is negative, the Schwarzian derivatives of all 
iterates are negative as well. 

Figure |5] shows the higher order iterates (of order 4, 8 and 16) in the same 
conditions as figure [T] We observe that the iterates become increasingly close 
together when their order increases, with smaller and smaller slopes at the 
fixed points corresponding to the 2-state regime. Moreover, they resemble more 
and more a square function, constant in various domains of the variable. This 
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Figure 8: Iterates for A = 0.95, 7 = 0.42 and C 0.8, of order 1, 2, 4, 8 and 
16. The convergence to the 2-state regime is visible. 



was expected: in the hmit of very large orders, whatever the variable is (i.e. 
whatever the initial conditions of the iteration are) one reaches a regime where 
only two values of the outgoing wave amplitude are possible; these values then 
remain stable, meaning that the action of more iterations will not change them 
anymore. So, one can read directly that the limit cycle is a 2-state on the shape 
of /^^, which has two values; it would for instance have 4 in the limit cycle 
was a 4-state regime for these values of the parameters. For the clarity of the 
figure, we have shown only iterates with orders that are powers of 2, but it is of 
course easy to plot all iterates. For a 2-state regime, even orders are sufficient 
to understand the essence of the phenomenon, since odd order iterates merely 

(2) (2) 

exchange the two fixed points and . 

In table 1, the existence of two different stable regimes for the same value 
of the parameters signals an inverse bifurcation; Figure [5] shows an example 
of such a situation. For 7 — 1.2, both the static and 2-state regimes are then 
stable, depending on the initial conditions. For the static regime, the curve f^^"^ 
coincides with the second diagonal y — —x, a case in which the fixed point is 
presumably stable (the stability becomes intuitive when one notices that the 
tangents of the higher order curves lie within the angle of the two diagonals). 
For the 2-state regime, the state of positive pressure value corresponds to a 
beating reed. 
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Figure 9: Iterates for A = 0.95, 7 = 1.5 and ( = 0.8, of order 1, 2, 3, 4, and 8. 
The curves of J*^®^ and f^-'-^^ are almost perfectly superimposed. Around x — 0, 
the convergence to the static regime appears to be very slow. On the contrary 
the convergence to the 2-state regime is rapid. 




Figure 10: Iterates for A = 0.95, 7 = 0.515 and C = 0.8, of order 1, 2, 4, 8 and 
16. 
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Finally FiglTU] shows another case of existence of two different regimes for 
the same value of the parameters. A 2-state regime can occur, as well as a 
4-state regimes can occur. It appears that the second one is more probable than 
the first one, when initial conditions are varied. 

5.2 Periodicity windows; intermittencies 

We now investigate some regimes occurring in a narrow range of excitation 
parameter 7. 

(i) We first examine a chaotic regime occurring just before a 6-state regime 
(period tripling) and the transition between the two regimes. Figure [TT] shows 
the iterated functions of order 1, 2, 6, and 12. The 6th iterated function crosses 
the first diagonal at the same points than the first and the second iterates only, 
which means that no 6-state regime is expected. By contrast, the 12th iterate 
cuts the diagonal at more points, but with a very high slope, indicating that 
the corresponding fixed points cannot be stable. This, combined with the fact 
that no convergence to a square function (constant by domains), such as /^^ 
in Figure [HI suggests an aperiodic behavior; the time dependent signal shown 
in FiglT2] looks indeed chaotic (nevertheless the flow always remains positive). 
The periodic/chaotic character of the signal can be distinguished by examining 
the time series, but a complementary method is the computation of an FFT. 
For the signal of FigfT^ the spectrum is more regular than the spectrum of a 
6-state periodic regime. Nevertheless the frequencies of the latter (the "normal" 
frequency /2 of the 2-state regime with the frequencies /2/3 and 2/2/8) remains 
visible in the spectrum of the first one, as it is often the case for signals cor- 
responding to very close values of the parameter. A consequence is that these 
frequencies clearly appear when listening the sound. 




Figure 11: Iterates for A = 0.95, 7 = 0.4445 and C = 0.8, of order 1, 2, 6 and 12. 
A convergence to an aperiodic regime is visible. The arrow indicates a region 
where f^^\x) is very close to the first diagonal, but does not yet cross it. 
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Figure 12: Iteration from n = for A = 0.95, ( = 0.8, 7 = 0.4445, = 0; the 
upper part shows the the pressure p, the lowest part the values of the flow u. 
The regime looks chaotic. 

Figure [T51 is similar to figure [TTl but with a slightly larger value of 7 (0.4469 
instead of 0.4445). In the region indicated by the arrow, one notices that the 6th 
iterated function now cuts the first diagonal. They are 12 points of intersection 
(plus 1 common point with the first iterate as well as two common points with 
the second iterate, all unstable); the slope of the tangent shows that 6 of them 
are stable, so that one obtains a 6-state, periodic, regime. The variations of 
higher order iterates, e.g. f^^'^\ remain very fast; the convergence to the limit 
cycle is then much slower than for Fig. |S1 except if the initial point is close to 
a limit point (e.g. that shown by an arrow: it turns out that the 12th iterated 
function is very close to the 6th one). As a consequence, the initial transient to 
the 6-state regime can be rather chaotic, as shown in Fig. [T31 but convergence to 
a periodic regime does occur later. This existence of periodic regimes above the 
threshold for chaos is called "periodicity windows" , which appears as a narrow 
whiter region in Fig. [5j A difference with the usual 2"-state regimes (when 7 
is below the chaotic range), for instance corresponding to Fig. [7l is that one 
obtains 2" intersections with the diagonal, stable or unstable; by contrast, for 
the 6-state regime, they are 6 stable and 6 unstable points. 

(ii) We now examine the transition between a 6-state regime and a 4-state 
regime through chaotic regimes or intermittency regimes. For 7 = 0.4544, a 
6-state regime is obtained. Fig. (TS] shows the iterates of order 1, 2, 4 and 6. 
The 4th and 6th iterates have common intersections with the first and second 
iterates, since both 4 and 6 are multiples of 2. The 6th iterate intersects the 
first diagonal at 12 other points, while the 4th cuts the diagonal at 4 points 
only. These 4 points are unstable, thus no 4-state regime can exist. On the 
contrary, for the 6th iterate, half of the 12 points are stable (i.e. with a small 
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Figure 13: Iterates for A = 0.95, 7 = 0.4469 and C = 0.8, of order 1, 2, 6 and 
12. A convergence to a 6-state regime is observed. The arrow indicates a region 
where /^^^(x) cuts the first diagonal. 

slope of the tangent line), so that one obtains a 6-state stable regime. 

What happens for a higher value of 7, namely 0.472 corresponding to a 4- 
state regime is shown in Fig. [161 with again the iterates of order 1, 2, 4, 6. 
The 4th iterate curve crosses the diagonal for the same number of points than 
previously, but the 4 points are now stable. The 6th order iterate does not 
intersect the diagonal, except at the common points with the two first iterates. 
Between the two preceding values of the parameter 7, both chaotic and inter- 
mittent regimes can exist. For 7 — 0.46623, Figure [17] shows intermittencies 
between a chaotic and a 6-state behaviors (upper curve), and Figure [TSj shows 
that the 6th iterate is tangent to the first diagonal in 6 points, so that the re- 
sulting permanent regime can be interpreted as a kind of "hesitation" between 
two behaviors. The 4 intersections of the 4th iterate remain unstable. 
The lower curve in Figure [T7] shows another, more visible, example of inter- 
mittencies, obtained with slightly different values of the parameters, between a 
chaotic regime and a 4-state one (actually it is a 8-state one, very close to a 
4-state regime). 
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Figure 14: Iteration from n for A = 0.95, C = 0.8, 7 = 0.4469, = -0.3347. 
The regime is periodic (6-state). 
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Figure 15: Iterates for A = 0.95, 7 = 0.4544 and C = 0.8, of order 1, 2, 4 and 
6. A convergence to a 6-state regime is observed. 
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Figure 16: Iterates for A = 0.95, 7 = 0.472 and C = 0.8, of order 1, 2, 4 and 6. 
A convergence to a 4-state regime is observed. 
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Figure 17: Iteration from n = for A = 0.95, C = 0.8, 7 = 0.46623, = 
(upper curve) : Intermittencics between chaos and a 6-statc! rcigime are observed. 
However the lower curve (for A = 1, ^ = 0.8, 7 = 0.467, = 0) shows a more 
clear situation of intermittencics between chaos and a 4-state regime. 
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Figure 18: Iterates for A = 0.95, 7 = 0.46623 and C = 0.8, of order 1, 2, 4 and 6, 
corresponding to intermittencies. The sixth iterate is tangent to the diagonal. 
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6 Conclusion 



The study of the iteration model of the clarinet should not be limited to the first 
iterate: higher order iterates give interesting information on possible regimes of 
oscillation. In the limit of very high orders, their shape gives a direct indication 
of the number of states involved in the limit regime, or of chaotic behavior. 
One can also predict an intermittent regime of the iterations, which takes place 
when an iterate is almost tangent to the first diagonal, so that the iterations are 
"trapped" for some time in a narrow channel. The phenomenon might be related 
to some kinds of multiphonic sounds produced by the instrument. It is true that 
this phenomenon takes place only in a rather narrow domain of parameters, but 
this is also the case of the period doubling cascade, which has been observed 
experimentally. One can therefore reasonably hope that the present calculations 
will be followed by experimental observations. 
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APPENDICES 



A Analytical iteration function 
A.l Derivation of the equations 

Our purpose is to obtain an analytical expression of the iteration function p+ = 
/(Pn-i)- From the basic model (Eqs. (|3] to [5j [71 [TT]) ) , the following quantities 
can be defined: 

X = 1 - Pn=l - Pti -Pn =1 - Pti + ^Pn-l'^ 

Pn = 5( P'n) be obtained from the knowledge of the function X{Y), given 
by the solving of: 

Y = X (beating reed, X > 1); (27) 

Y = a: + C(l -^)\/^ (non-beating (28) 

reed, positive fiow, < A < 1); 

Y = A - C(l - A)V^ (non-beating (29) 

reed, negative flow, A < 0). 

For the non-beating reed case, the study of function Y{X) leads to a direct 
analytical solution, as explained below, at least if ^ < 1 (otherwise it is a multi- 
valued function). 
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Finally, with the notation x — p^-i and f{x) — p^, if y{X) is the Heaviside 
function, the iteration function is obtained, as: 



f{x) = 7-X(r) + Aa;, withy = 7 + 2Aa; and 
Y{X) = X + Csi9n{X)y{l- X){1- X)y/\X\. 



(30) 
(31) 



A. 2 Non-beating reed, positive flow (0< Y < 1) 

For this case, both X and Y are positive and smaller than unity, because C < 1- 
Writing Z = \fX, Eq. (gS]) is written as: 



Gx{Z) = Y, where Gi{Z) = -(Z 



Z^---l 

C 



(32) 



The study of function Gi{Z) shows that it is monotonously increasing from to 
1 when Z increases from to 1. Therefore the equation Gi{Z) — Y has a unique 
solution when < F < 1. With this condition, it appears that the equation has 
three real solutions, and that the interesting solution (located between and 
1) is the intermediate one. As a conclusion, it is possible to use the classical 
formula for the solution of the cubic equation: 



X 



Z 



-T] sm 



1 



tfj ; 



V' - 
9 



1 

3C 



!)• 



A. 3 Non-beating reed, negative flow (Y< 0) 

For this case, both X and Y are negative. Writing Z — \/ 
written as follows: 



-X, Eq. 



G2{Z) = y, where G2{Z) = -C^Z 



(33) 



The study of the function Gi (Z) shows that it is monotonously decreasing from 
when Z increases from 0. Therefore the equation Gi{Z^ = Y has a unique 
real, positive solution when Y <Q. The two other solutions are either real and 
negative or complex conjugate, with a negative real part, because the sum of 
the three solutions is negative (— 1/C). As a conclusion, the solution can be 
written by using the following formulae: 

If the discriminant is positive 

discr ~ + > 0, where 




[3-^] 



V^ + M 
27C 



9 

Z = Si W7'^ si = 

si 3C 



Vdi 



1/3 
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If the discriminant is negative 

discr = q'^ + < 



Z = -77' cos 



— arccos 



+ M 

' Of 



1 

3C 



B Negative flow limit 

The condition of existence of negative flow is given by Xi > /min- This is equiv- 
alent to the condition on the antecedents, x'^ < /max, where x'^ is the larger 
antecedent of Xi, such as x[ > x^ax, because f{x) is decreasing for all x > Xmax 
(see Fig. [5]). Therefore the volume flow is negative at time n + 1. 

In order to determine the limit value 7^, the following equations are to be 
used: 



Y = j + 2Xx'^^j{l + X) + X2A^. 



(34) 
(35) 



7 being positive (a reasonable hypothesis for the normal playing), the unknown 
X needs to be larger than the quantity XA,^. Eliminating 7 in the above equa- 
tions implies the following equation, with X > XA(^: 



Y{X) ^X = 



{X-1)X + 2XA1; 



l + X 



or 



HiX) = S 



(36) 



with: 



H{X) = {1 + X)[Y{X)-X] + {1-X)X (37) 
for X > XA^ ; 6^2XA^. 

An example of function H{X) is shown in Fig. [121 It appears that no solutions 
exist if H{1) > 6 and two solutions exist if H{1) < S, i.e. if inequation ([M)) holds. 
The two solutions can be obtained analytically. However, for sake of simplicity, 
we give the exact solution for the larger one, 7^2, and an approximation for the 
smaller one, 7^1, obtained at the first order in s — 1 — X: 

7i2 = 

7^1 - 



2X 
with e = 



[1 - AAc - e] , 
A - 1 + 2AA, 



C 



(A + 1)C + A- f 



(38) 
(39) 
(40) 
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Figure 19: Function H{X) given by Eq. ([57]) and constant line 6 — 2Ayl^. Two 
solutions X > XAi^ exits for this case (A — 0.9, ( — 0.9), because condition 
is satisfied. 

This error is found to be less than 1% in comparison with the exact value. 
Condition ((24)) can be shown to be necessary and sufficient. We do not give the 

entire proof, but it can be shown that another necessary condition for having 
two solutions is H'{1^) < 0, or ^(A + 1) + A — 1 > 0, but it is implied by 
condition (OI)) . 

Fig. |4] shows that the first negative fiow threshold 7^1 is very close to the 
threshold 7b, and slightly smaller. For a given A, the limit value of C such as 
A > 1/(1 + 2Aq) corresponds to the equality between the beating reed threshold 
and the negative fiow one. For a given negative fiow is possible above a certain 
value of A. For rather strong losses, if A < 0.84, no negative fiow can occur. For 
a cylindrical resonator, this implies that ai > 0.085. 
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